One-magnon Raman scattering in L^CuO/i: the origin of the field-induced mode 
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We investigate the one-magnon Raman scattering in the layered antiferromagnetic La2Cu04 com- 
pound. We find that the Raman signal is composed by two one-magnon peaks: one in the Bi g 
channel, corresponding to the Dzyaloshinskii-Moryia (DM) mode, and another in the B-i g chan- 
nel, corresponding to the XY mode. Furthermore, we show that a peak corresponding to the XY 
mode can be induced in the planar (RR) geometry when a magnetic field is applied along the easy 
axis for the sublattice magnetization. The appearance of such field-induced mode (FIM) signals 
the existence of a new magnetic state above the Neel temperature Tn, where the direction of the 
weak-ferromagnetic moment (WFM) lies within the Cu02 planes. 
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The possibility of inelastic light scattering by one- and 
two-magnon excitations in magnetic insulators was ac- 
knowledged long ago by Elliot and Loudon pj and the 
microscopic theory for suchprocess was described in de- 
tail by Fleury and Loudon |2J]. One of the possible mech- 
anisms for the magneto-optical scattering is an indirect 
electric-dipole (ED) coupling via the spin-orbit interac- 
tion, and such mechanism has been used to determine 
the spectrum of magnetic excitations in many different 
condensed matter systems like fluorides, XF2, where X 
is Mn 2+ , Fe 2+ , or Co 2+ Q, inorganic spin-Pierls com- 
pounds, CuGe03 0|, and the parent compounds of the 
high-temperature superconductors, such as La2Cu04 0. 

In this communication we report on a theoretical study 
of the Raman spectrum in a two-dimensional (2D) quan- 
tum Heisenberg antiferromagnet (QHAF) with DM and 
XY interactions, which is a model for the magnetism of a 
single layer in La2Cu04. As it was recently shown in 
although small, such anisotropic DM and XY terms are 
very important in order to explain the unusual magnetic- 
susceptibility anisotropies observed in La2Cu04 [f|. It 
was found, in particular, that, below the Neel ordering 
temperature, Tjv, the two transverse magnon modes are 
gaped, with their respective gaps being uniquely deter- 
mined by the strength of the DM and XY interaction 
terms, in agreement with neutron scattering experiments 
Q. However, more recent Raman-spectroscopy experi- 
ments reported on the presence of only one of the above 
two magnon modes, the DM mode, in the Raman spec- 
trum for the planar (RL) geometry (i?=right-rotating 
and L=left-rotating) at zero magnetic field 0. Most 
surprisingly, it was also found that a second mode can 
be induced in the forbidden (RR) geometry, for a mag- 
netic field applied along the easy axis for the staggered 
magnetization, and only for such configuration. As we 
now clarify, the appearance of such FIM is directly asso- 
ciated to a continuous rotation of the spin quantization 



basis, which modifies the selection rules set by the ED 
coupling on the Raman scattering cross-section, thus al- 
lowing for the observation of the, previously selected out, 
XY magnon mode in the spectrum. 

The Hamiltonian representing the interaction of light 
with magnons can be written quite generally as || 



Bed = ^E£ X (r)E/, 



(1) 



where Eg and Ej are the electric fields of the scattered 
and incident radiation, respectively (a T is the transposed 
of the a vector) and x( r ) is the spin dependent suscep- 
tibility tensor. We can expand \( r ) m powers of the 
spin-operators, S, as 
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where [i,v = x,y,z label the spin components. The low- 
est order term Xg^(r) is just the susceptibility in the ab- 
sence of any magnetic excitation (it corresponds to elas- 
tic scattering), and it will be neglected in what follows. 
The second and third terms can give rise to one-magnon 
excitations because they can be written as S ± (r) and 
S z (y)S ± (y) respectively. The intensity of the scattering, 
as well as the selection rules, will be determined by the 
structure of the complex tensors K and G. For a square- 
lattice antiferromagnet the linear term on the spins re- 
duces to (for Stokes scattering only) 



E t E i)S~ A 



~( E s E i - E S E l) S i_ 



(3) 



where the sum runs over the two sub-lattices A and B, 
the Ko coupling constant is the same for the two sub- 
lattices, and z is the direction of the spin easy-axis jjj. 
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Analogously, the quadratic term on the spins reduces to 
H S ED = G Q J2 K E s E t + E s E ?)Oi A 

iA,lB 

+ (E z s Ey +E s Ef)0+ B ], (4) 

where Go is a coupling constant and Of — SfS\ +S i Sf 
||. We see that while the K tensor is purely imaginary 
and totally anti-symmetric, the G tensor is purely real 
and totally symmetric, as required by the general sym- 
metry properties of the electric susceptibility of a mag- 
netic material, also known as Onsager relations Q. Fur- 
thermore, the relative values of the coupling coefficients, 
Kq and Go, can be deduced from the measurements of 
magneto-optical effects like: magnetic circular birefrin- 
gence (to determine Kq) and magnetic linear birefrin- 
gence (to determine Go) ||. 

The total ED Hamiltonian H E d = Hed + ^ed can 
be finally written in terms of the x,y,z components of 
the sublattice magnetization, Mi = (Si A — Sj B )/2, as 

ff Br > = £ {E£ X *E, Mf + V T sX v ViM?}, (5) 
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where we introduced the matrices 
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with d = GqS + iKq. Here we made the usual mean- 
field assumption (Sf ) = — (S? B ) = —S and we dropped 
terms of the type Sf* + Sf* since these will give rise 
to negligible contribution for the k = scattering in the 
long- wavelength limit. In deriving Eq. © we assumed 
that the spins order along the z direction. In the more 
general case, the spin easy-axis is along a different direc- 
tion z with respect to a given (xyz) reference system (for 
example, the one attached to the unit cell of the crystal). 
In this case, one should rewrite the ED Hamiltonian as 

H ED = { E sf E i % + ^sX v ^i M?} , (6) 

i 

where X x ' v = R T X X ' V R, M = i?M, and R is the matrix 
of the rotation from (xyz) to (xyz). 

The above Hamiltonian (§J is to be treated as a pertur- 
bation of the following single-layer 5 = 1/2 Hamiltonian 
for the La2CuC>4 system 



■E 



• (Si x Sj) + 2__^ • r ij ■ Sj, 

(7) 

where J is the in-plane antiferromagnetic super-exchange 
between the spins, Si, of the neighboring Cu 2+ ions, 
and Dij and T y are, respectively, the DM and XY 
anisotropic interaction terms of the low temperature or- 
thorhombic (LTO) phase of La 2 Cu0 4 OH- The long- 
wavelength effective theory for the Hamiltonian {JJ was 
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FIG. 1: Left: LTO coordinate system with respect to a single 
Cu02 layer. Right: direction of the staggered magnetization, 
(n), of the DM vector, D+, and of the WFM, (L) oc (n)xD+, 
for B = 0. In the scattering geometry of Q the electric field 
vector of the light, E/,f, is always parallel to the ab plane. 



derived in Q. As usual, one decomposes the Cu 2+ spins 
in their staggered, n, and uniform, L, components as: 
Si(r)/S = e iQ ' x *n(xi,T) +L(x i ,r), where Q = (vr,7r), 
and one then integrates out L to obtain 

S = — !— [ dr [ d 2 x{(<9 r n + iB x n) 2 +c 2 (Vn) 2 
2.9oCo Jo J 

+ (D + -n) 2 +r c (n c ) 2 + 2B-(D + xn)}, (8) 

where n(x, t), the continuum analogous of Mi, is re- 
quired to satisfy n 2 = 1. Here go is the bare coupling 
constant, related to the spin- wave velocity, cq, and renor- 
malized stiffness, p s , through p s = Co(l/Ngo — A/4ir) jjj, 
A is a cutoff for momentum integrals (we set the lattice 
spacing a = 1), N = 3 is the number of spin compo- 
nents, D + = -D+x a , and T c > 0. From now one we will 
use the LTO (abc) coordinate system of Fig.QJ where x a , 
Xb, and x c are the LTO unit vectors. From Eq. JHJ it 
follows that the two transverse staggered modes n a and 
n° have their gaps given by A a = D + = 2.5 meV and 
A c = VTc = 5 meV, respectively, so that Xb is the easy- 
axis for the staggered magnetization. Since the uniform 
magnetization (L) = (1/2 J) ((n) x D+) @], we find that, 
in the ordered AF phase (for B = 0), a non-zero weak- 
ferromagnetic moment (WFM) (L) is present, directed 
along x c , see Fig. 1 and Fig.E^a). 

We can now discuss the experiments of 0] in the ab- 
sence of magnetic field, B = 0. First, let us observe that 
since the sublattice magnetization, (n), is oriented along 
the LTO Xb direction, the Hamiltonian © reads in the 
continuum 

H ED = J d 2 x {E^ X a E/ n a + e£ X c E 7 n c ) 

= y"d 2 x{n a n a +n c n c }, (9) 

where we have introduced the polarization projectors 
II a , c = X a ' c Ej with 



(10) 
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in the LTO abc coordinate system. The structure of the 
above matrices implies imediately that the n a (DM) and 
n c (XY) modes should be observed, respectively, in the 
Big and B 3g channels of the non-magnetic Dmab or- 
thorhombic group of the LTO phase of La2Cu04. 

The Brillouin zone center one-magnon Raman inten- 
sity can now be calculated from Fermi's golden rule using 
the Hamiltonian JjjJ as a perturbation, and we obtain, for 
Stokes scattering, 

= [n B (w) + 1] ||n Q | 2 A(o,^) + |n c | 2 A(o,^)} , 

(ii) 

where 71b(o>) = (e' 3 " — 1) 1 is the Bose function and 
A a , c (0, oj) is the q = spectral function of the transverse 
modes, which is peaked at the mass value A QjC . 

The scattering geometry used in is the so called 
backs cattering geometry. In this setting, the direction of 
the propagating wave-vectors of both the incoming and 
outgoing radiations are chosen to be along x c , in such a 
way that the electric field of the light is always parallel 
to the Cu02 or ab plane, see Fig. ^ Let us denote by 
(einZout) the polarization configurations. Two polariza- 
tion configurations were used in jj|: i) light linearly polar- 
ized with e-in = x Q and e Q ut = *-b (Bi g geometry); and ii) 
light Right/Left circularly polarized, e R — -7= (x a — ix b ), 

and e L = ^(x a + ix b ), combined in the (RL) or (RR) 
geometries. It is worth mentioning at this point that 
the (RL) and (RR) polarization configurations probe, 
respectively, the anti-symmetric (imaginary) and sym- 
metric (real) parts of the Raman tensor \ a m HIM , i.e., 
U R L = 2iK and U RR = 2iG S. Moreover, because the 
electric-field of the light is always parallel to the ab plane, 
n c = for any polarization configuration. As a result, 
we find that only the spectral function of the n a mode, 
A a (0, u>), which is peaked at the DM gap, A a , contributes 
to the Raman intensity 11 111 in the Bi g channel, in agree- 
ment with the experiments of 0| for both the B\ g and 
(RL) geometries. We observe, furthermore, that within 
the ideal microscopic theory of Fleury and Loudon Q, 
one should always have Gq = for magnetic systems in 
which the ground state has zero or quenched orbital an- 
gular momentum. When this happens, also H RR = 0, 
in which case we could refer to the (RR) geometry as 
forbidden. 

Let us discuss the effects of a magnetic field B applied 
along x b , which is the configuration where the FIM is 
observed in Ref. Qj]. The first term in Eq. © is re- 
sponsible for the softening of the transverse gaps, by an 
amount B 2 , A 2 — > A 2 c — B 2 , as observed in the mea- 
surements of Ref . 0| . However, this softening is quanti- 
tatively very small for small fields and will be neglected 
in what follows. Conversely, the most remarkable effect 
of the magnetic field comes from the last term in Eq. 
JHJl, which can be written as (l/goco) J Q dr J d 2 x h n c , 
where h = |B x D+|. As we can see, such term gener- 
ates, for B || Xf,, an effective staggered field h coupled to 
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FIG. 2: Orientation of the staggered magnetization (n) and 
of the WFM (L) for: a) zero field B = 0; b) B < B c and 
T < Tjv (8 and (n) are functions of T); and c) B < B c and 
T > Tm. Observe that c) is also realized at any T for B > B c . 



n c , leading to a rotation of the staggered-magnetization 
order parameter with respect to the zero-field case, see 
Fig. EJb). In fact, if we write (n) = (0,<Ti,,a c ) for the 
order parameter we obtain at mean-field level 



l-a 2 -NI ± (^^), a c 



-h/r c 



(12) 



for T < Tn and 



a b = 0, <7 c = -/i/(r c + r 2 ), 1 = <t 2 c +NI ± (£), (13) 

for T > Tn, where, following the notation of Ref. 
[si, £ is the correlation length and I± = (1/2) (la + 
I c ) is the integral of the transverse fluctuations |l2| 
I a = ( 5 r)(2^c)log{(sinh( C A/2T))/(sinh(A Q /2r))}. 
The temperature dependence of the components of the 
(n), as well as of the angle 6, 



tan 9 



(14) 



that defines the rotation of (n) from the Xf, direction, are 
reported in Fig. Observe that for h 7^ we conclude 
that the uniform magnetization 



£± 

2J 



(|(T b |x c + |(T c |x b ) . 



(15) 



is also rotated from its original orientation, perpendicular 
to the Cu02 planes, by the same angle 9, see Fig. |2fb). 
Notice, furthermore, that such rotation has the x a axis 
as the symmetry axis, in such a way that both (n) and 
(L) are always confined to the bc plane 0. Moreover, 
the angle of rotation varies continuously from its T = 
value until 9 = — tt/2, see inset of Fig.03 because cr b — > 
continuously as T — > Tn(B). This is a very important 
finding because it shows clearly that no two-step spin-flop 
transition occurs, in agreement with magnetoresistence 
experiments ^lj. Finally, as the field B increases the 
T = value of cr b decreases, until a critical field B c above 
which the order parameter is always oriented along the 
x c direction. As a consequence, either at T > 2jv(B), 
for B < B c , or at arbitrary temperature and B > B c 
we find a new magnetic state with cr b = and a c 7^ 0, 
where the WFM is confined to the ab plane and ordered 
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along the Xf, LTO direction, as it was suggested by the 
experiments of Q • Within our mean-field calculation the 
critical field B c at T = can be estimated from Eq. 1(12)1 
as B c = yjl - NgA/A-K(T c /D + ) = 40 T. 




-0.1 



0.1 0.2 0.3 0.4 0.5 
T/J 



0.1 



0.2 0.3 
T/J 



0.4 



0.5 



FIG. 3: T dependence of a b and o~ c at H = 10T, from Eqs. 
112t and 1131 . Inset: T dependence of the angle 9, from Eq. 

We are now ready to explain how the above rotation of 
the spin easy axis, induced by the magnetic field, mod- 
ifies the Raman spectra. The matrix R describing the 
rotation by the angle 9 around the x Q direction modifies 
the Raman matrices iflOjl to 



this case, (n) = ObXb, and (L) = (1/2 J)D + ai,x c , see Fig. 
I^a). The case of B || x c was already discussed in 0|. Al- 
though in this third case h ^ 0, it is actually a source for 
the n b component of the staggered magnetization, and 
we again have a c = 0. Thus, for either B || x a or B || x c , 
where 9 — 0, we will have II C = and only the DM mode 
will be present in the one-magnon Raman spectrum. 

In conclusion, we have found that the appearance of 
the FIM for B || x;, is a consequence of both a continuous 
rotation of the spin-quantization basis and of a nonzero 
symmetric part for the \ c Raman tensor, Go ^ 0, signal- 
ing a deviation from the ideal ED scattering mechanism 
of Fleury and Loudon Q. Furthermore, we have found a 
new magnetic state, well above TV, in which the WFM is 
confined to the ab plane. The absolute intensities of the 
one-magnon Raman peaks in the two circular polariza- 
tions, (RL) and (RR), are given in terms of K and Go, 
which have to be determined from first principles. More- 
over, to obtain the temperature and field dependencies of 
the relative peak intensity of the FIM a more precise cal- 
culation of the T and B dependencies of the anisotropy 
gaps and of the spin damping would be required. 

The authors have benefitted from invaluable discus- 
sions with Y. Ando, G. Blumberg, S. Caprara, A. Gozar, 
M. Grilli, A. N. Lavrov, and J. Lorenzana. 
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The new transverse spin modes correspond to h a = n a , 
since the rotation leaves the x a direction untouched, and 
n° = n c cos 9 + n b sin 9. Observe that the main effect 
of the rotation is that now II C cx E^x c E/ ^ 0, so that 
also the n c (XY) mode can be observed. More specif- 
ically, from Eqs. © and ((Tfil) it follows that in the 
(RR) geometry used in Ref. jj| U RR = 2zGo cos 9 and 
H RR = Go sin 29 cos 9. As a consequence, now also the 
spectral function of the n c mode, A c (0,lu), contributes 
to the intensity 11 111 with a peak at the energy A c , which 
corresponds to the XY gap. We thus conclude that the 
FIM mode observed in Ref. 0] is nothing less than the 
XY mode, which can only be seen for B || Xb and nonzero 
Go. We should remark at this point that the above effect 
occurs only for B || Xf,, thus justifying the appearance of 
the FIM only for that orientation of the magnetic field. 
In fact, for B II x a , we have h = and <r c — 0. In 
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